Ramsey interferometry using the Zeeman sublevels in a spin-2 Bose gas 
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We perform atom interferometry using the Zeeman sublevels of a spin-2 Bose-Einstein conden- 
sate of 87 Rb. The observed fringes are strongly peaked, and fringe repetition rates higher than the 
fundamental Ramsey frequency are found in agreement with a simple theory based on spin rota- 
tions. With a suitable choice of initial states, the interferometer could function as a useful tool for 
magnetometry and studies of spinor dynamics in general. 



The study of spinor properties of matter was extended 
to the field of dilute Bose-gases after the achievement 
of optical confinement of a Bose-Einstein condensate 
(BEC) p. This technique gave rise to a number of in- 
triguing studies at the intersection of fields as diverse 
as magnetometry, symmetry breaking and pattern for- 
mation [2]-[4]. More recently, spinor [5] and effective 
spinor [6 condensates (i.e. two mode condensates which 
can formally be treated as spin 1/2 systems) along with 
atom-atom interactions have provided the tools to cre- 
ate "non-linear" atom interferometers, that is, interfer- 
ometers where noise in the quantity of interest can be 
suppressed by using spin-squeezed states [TJ [8] , leading 
to measurement sensitivity in excess of that classically 
permitted. 

Spin rotations for spin > 1/2 have been made use 
of for decades in molecular experiments [9]. Moreover, 
in terms of interferometry experiments, sharply peaked 
fringes were observed in an atom interferometer using 
magnetic sublevels in [10 . In that study the source of 
the fringe pattern was explicitly linked to multi-path 
interference. Additionally, Ramsey interferometry in a 
spin-1 atomic gas has been used in sensitive magnetom- 
etry experiments, with the aim of detecting a permanent 
electric dipole moment [11] . Of relevance to our particu- 
lar experiment, we note that rotation and more general 
manipulation of spin states have been used as tools in 
various studies of spin-2 spinor BECs (see for example 

0CE2CL3]). 

Here, we introduce a spinor Bose-gas interferometer 
using a spin-2 condensate of 87 Rb, in line with recent 
theoretical proposals p~4j [15] . We observe predicted non- 
trivial fringe structure, with sharp fringes being a partic- 
ular feature of interest [14], and, as we will show, fringe 
oscillations at multiples of the fundamental period being 
another. Parallels with proposals for spin echo experi- 
ments in spinor BECs can also be made [IB] . We note 
that our work shares common tools and certain similar- 
ities with the experiments mentioned in the above para- 
graph, but our principle aim is to demonstrate the pos- 
sibilities of Ramsey type interferometry in a spin-2 BEC 
setting. We also note that while non-destructive imag- 
ing techniques can extract precise information about the 
atomic precession frequency [2], performing an interfer- 
ometer sequence followed by a full Stern- Gerlach (SG) 
measurement as we do here can provide more informa- 




time 



FIG. 1. Schematic diagram of the experiment, (a) Depicts 
the BEC confined in a crossed-dipole trap where the inter- 
ferometer sequence takes place. The arrow below the BEC 
indicates its movement under gravity after it is released from 
the trap, and its subsequent splitting into five components (la- 
belled vriF — —2,. ..,2) due to an applied gradient magnetic 
field, (b) shows the interferometer sequence in more detail. 
Grey arrows depict the atomic spin, while the RF field is indi- 
cated by a red arrow. Note that the axis of the bias magnetic 
field is shown by a thin green line labeled "B" in both (a) and 

(b). 



tion about the system. 

The concept of the experiment is illustrated in Fig. [I] 
The five magnetic sub-level components associated with 
a spin-2 gas can be separately measured, and the pop- 
ulation of each component depends on coherent rota- 
tions applied by a radio frequency (RF) field and also 
Larmor precession. In our experiment, a BEC of 87 Rb 
is created in a magnetic trap using standard RF evap- 
oration methods. After loading into a crossed optical 
dipole trap, typically 4xl0 5 atoms remain, and the trap- 
ping potential has axial and radial frequencies of order 
~ 33Hz and ~ 100Hz respectively. The spin orienta- 
tion of the optically trapped condensate is preserved by 
application of a 20G bias magnetic field (B-field) and 
the optically trapped BECs internal state is given by 
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FIG. 2. Coherent spinor rotation of the BEC. (a) - (e) show 
the behavior of P — |^/| 2 for the ttif components +2,-+,0,- 
1,-2 respectively. In each case squares show the experimental 
data and dashed lines show the theory of Eq. [I] (f) and (g) 
show the effect of changing the RF frequency on column den- 
sities p for a nominal tt/2 RF pulse of 29/xs in width for the 
near resonant and off-resonant cases respectively. Finally, (h) 
shows the component- wise distribution of P = |^| 2 for a de- 
tailed scan over the resonant frequency. The populations of 
the wif = +2, +1,0,-1 and —2 states are shown by solid, 
dashed, dash dotted, fine dotted and large dotted lines re- 
spectively. 



\F = 2 1 mp = +2). The mp components may be sepa- 
rated using a SG technique as shown in Fig. [TJa). After 
a 200ms hold time in the optical trap, we ramp the bias 
field down to 300mG over 10ms, and apply the interfer- 
ometer sequence depicted in Fig. [TJb). The 300mG field 
is applied by separate coils, the current supply for which 
is a commercial laser diode supply with current ripple on 
the order of 10/iA. We found that a stable bias field was 
very important for achieving repeatable spin rotations. 

During the interferometer experiment, an RF tt/2 
pulse rotates the atomic spin until it is perpendicular 
to the bias B-field. Subsequently, Larmor precession is 
allowed to occur for a time T. We calculated the Larmor 
frequency to be fo = gpjiBB/h = 210kHz at our esti- 
mated bias B-field of 300mG, where gp is the g- factor, 
and hb is the Bohr magneton. Finally a second RF tt/2 
pulse resonantly rotates the spin again, and the popu- 
lation of each spin component is read out using the SG 
technique followed by absorption imaging of the conden- 
sate. Further details of the experimental setup may be 
found in Refs. [13 [18]. 

Before presenting our experimental results, we offer 
a phenomenological theory of the experiment described 
above. Our experiment realizes a Ramsey- type interfer- 
ometer over the five Zeeman sublevels of our 87 Rb con- 
densate. The usual treatment of such systems assumes a 
two level system and indeed, the theoretical treatments 
of spinor BEC interferometers so far transform into an 
effective two level system using the Major ana represen- 
tation [151120]. 



Here, we adopt a different approach, and derive a for- 
mula for the interferometer fringe (i.e. population as a 
function of RF wave frequency) in each sub-level. We 
rely on the following two facts: (i) An RF pulse with 
a frequency near to the Larmor frequency of the atoms 
induces a coherent rotation of the atoms about an axis 
perpendicular to the B-field [21-23 and (ii) the rotation 
effected by the RF pulse drops off as a smooth function 
of the RF detuning from the Larmor resonant frequency. 

We tested both of these assumptions experimentally 
and the results are shown in Fig. [2] In Fig. [2|a-e), mea- 
sured populations (squares) in mp sublevels +2, . . . , — 2 
respectively are compared with theoretical values. The 
theoretical values may be obtained by recalling that for 
non-trivial rotations j3 (i.e. rotations which couple rnp 
sublevels) the Wigner rotation matrix using the Euler 
convention is, D 2 (a = 0, /?, 7 = 0) = d 2 ((3) where d 2 (fl) is 
the so-called reduced rotation matrix [24] . Applying this 
operation to our initial state (all atoms in the mp = +2 
sublevel) * tf = [^ + 2,^+i,^o,^-i^-2] T = [1,0,0,0,0] T 
and squaring gives the measured probability distribution 
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For comparison with experiment we take (3 = ut, where 
t is the RF pulse width and u is the angular rotation 
frequency of the spin induced by the RF pulse. As seen in 
Fig. [2ja-e), the experimental values agree very well with 
theory for a measured value of u = 2tt x 8.8kHz. We 
can therefore be confident that the RF pulses perform 
coherent rotation of the spin as expected, and that the 
second order Zeeman effect is negligible on the time scale 
considered here. In that case, the five separate fringe 
patterns produced (one for each Zeeman sub- level mp) 
are non-sinusoidal, and their shape is predicted by the 
properties of the Wigner rotation matrix D 2 [24 . The 
peak structure can be explained by invoking the spin-2 
rotation properties of the BEC, rather than an explicit 
multi-path interferometer formulation. 

Regarding our second assumption, Figs. |2jf-h) show 
the effect of a nominal tt/2 pulse as the frequency is swept 
over the resonant Larmor frequency /q. The rotational 
response of the spin induced by the pulse is seen to tend 
to zero as the size of the detuning of the RF frequency 
/ from fo increases. The width A of the frequency re- 
sponse is seen to be A ^ 25kHz. By analogy with the 
standard Ramsey interferometer, we model the frequency 
dependence by a sine function. In particular, 
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The frequency dependence of the interferometer tt/2 
pulses can now be modeled by the operator d 2 [(7r/2)rj(f)] . 
Applying this operator to the initial state models the 
application of the interferometer's first tt/2 pulse. 
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FIG. 3. Measured interference fringes. Interference fringes are shown for a pulse separation T = 285/xs. (a) shows raw column 
densities p from absorption images of the BEC. The population is seen to oscillate between the = ±2 states as a function of 
frequency with the amplitude flowing through the tuf = +1,0, and —1 states, (b-f) show relative population measurements P 
for ttif — +2, . . . , — 2 respectively, with discrete points showing the experimental results (each being an average over the results 
at three neighboring frequencies) , error bars showing the standard deviation over these three points and solid lines showing the 
fitted theory from Eq.[3]in each case. The fitted parameters are fo = 195kHz and T = 290/xs and the peaks have a phase offset 
of 0.14 radians. 



To model the entire interferometer, we also need 
to model the phase accumulation due to Larmor pre- 
cession between pulses and the phase difference be- 
tween the RF pulses after the interrogation time T has 
elapsed. Both of these effects can be modeled by ro- 
tations about the B-field axis by applying the rotation 
matrix D 2 ($(f,T),f3 = 0,7 = 0), where $(/,T) = 
2n(f + fo)T is the combined phase accumulation due to 



Larmor precession and the time between pulses. Since 
f) = 0, the rotation matrix is diagonal with matrix el- 
ements D^ TlF mF = exp(— irriF&). Finally, the rotation 
operator d 2 [(it/2) rj(f)] is applied again to give the final 
output state. To evaluate the interference fringe, we ap- 
ply the Born rule to the complex amplitudes for each 
tuf component to find the relative population |\I//| 2 = 
|d 2 [(7r/2)7 ? (/)]^ 2 ($,0,0)(i 2 [(7r/2)7 ? (/)]^ i | 2 . Performing 
the matrix multiplication gives 
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where we have left the frequency dependence of and 
r] out of the notation for brevity. For comparison with 
experiments, we allow an additional free phase parameter 
(j> in the expression for <£. This phase depends linearly 
on T alone and its effect is to shift the position of the 
peaks without changing their spacing. 



We now move on to our experimental interferometry 
results. To perform interferometry on the BEC, we fixed 
a value T for the time between the RF pulses (the in- 
terferometer interrogation time) and observed the rela- 
tive population of the Zeeman sublevels as the pulse fre- 
quency / was varied between 175 and 210kHz. Fig. [3ja) 
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FIG. 4. The effect of interrogation time and input state. 
Fringes over approximately one period showing the effect of 
different interrogation time T. (a) shows the populations for a 
nominal pulse separation of T = 152/is, with mf as indicated 
on the vertical axis label, (b) shows data for a nominal pulse 
separation of T = 530/xs. Note that the horizontal scales in 
each plot are the same for ease of comparison. As in Fig. [3j 
discrete points show experimental measurements (this time 
without averaging) and lines show a fit of Eq. [3] to the data. 
The fitted parameters are fo — 207kHz and Tfi t = 143/is for 
(a) and fo — 192kHz and Tfit = 574/xs for (b) and the fringes 
were found to have phase offsets of -0.97 and 2 radians for 
(a) and (b) respectively, (c) and (d) show the effect of much 
longer interrogation times where phase changes due to effects 
other than rotation cannot be neglected, (c) shows the col- 
umn density after a single tt/2 pulse and 5ms hold time in 
the trap, (d) shows the column density for the same situa- 
tion as (c) but with a final tt/2 pulse applied at the end of 
the sequence. Note that the essentially uniform distribution 
of amplitude among states cannot be arrived at by simple 
rotations. Lastly, Figs, (e-g) show how the interferometer 
can be improved by using a different initial state, (e) and (f) 
show fringes for the ttif — +2 and tuf = states respec- 
tively, in the usual case where the initially populated state is 
ttif — +2. (g) shows the tuf = fringe for the case where 
the initial state is \rriF = +1)- In this case, we see that 4 
times as many fringes per period are expected compared with 
the fundamental periodicity of the interferometer. The gray 
shaded region shows one period (i.e. 2tv phase change) in each 



shows normalized column densities calculated from ab- 
sorption images of the BEC for a range of frequencies. 
From this essentially raw data, the oscillatory behavior 
is clearly visible, but the detailed structure is not read- 
ily apparent. To see the fringe structure, we calculated 
the relative population in each substate (shown as dis- 
crete points in Figs. [3^b-f)). We then fitted Eq. [3] to the 
experimental result for mp = +2 ( which has the best 
signal to noise ratio) and used the obtained parameters 
in Eq. [3] for all mp values. The fits are shown by solid 
lines in Figs. [3^b-f). The fitted parameters agree with 
the experimentally measured parameters to within 3% 
with a small extra phase shift of (j) = 0.14 radians being 



necessary to account for the exact peak position. 

The experimentally measured fringe pattern is non- 
sinusoidal and in good qualitative agreement with the 
predictions of Eq. [3j It is interesting to note that the 
mp = +1,0 and —1 component fringes contain oscilla- 
tions with a smaller period than the fundamental seen in 
mp = ±2. In particular, the mp = fringe oscillates 
with essentially half the period of the principle oscilla- 
tion. This kind of behavior is also found in multi-beam 
optical interferometers [25] and is a logical consequence 
of the fact that each peak in the population of mp = ±2 
is defined by two zero-crossings (i.e. as the amplitude 
crosses through the intermediate mp states). 

We estimated the sensitivity of the interferometer by 
finding the standard deviation over three nearest neigh- 
bor frequencies for every second point in our frequency 
scan. Then, for each peak in the highest visibility 
mp = +2 fringe pattern, we found the smallest distin- 
guishable phase difference within the experimental error. 
In this way, we calculated the average phase sensitivity 
to be about 0.6 radians or 10% of a fringe. Given our 
relatively short interrogation time, this result is compa- 
rable with other Ramsey atom-interferometers operating 
on coherent states [26]. We would like to note, however, 
that utilizing all five fringes of our spin-2 system presents 
an advantage relative to spin 1/2 or spin 1 systems be- 
cause at any given frequency five fringes are available. 
By choosing the fringe with the steepest rate of change 
at the desired frequency, sensitivity can be maximized. 
Additionally, we note that engineering the initial distri- 
bution of mp states can produce high visibilities in states 
other than mp — ±2. 

We also tested the interferometer for shorter and longer 
interrogation times, as seen in Figs. |4^a) and (b). The 
fringes show typical Ramsey-type behavior with the pe- 
riod becoming longer or shorter as the interrogation time 
was decreased or increased respectively. Details are given 
in the caption of Fig. [4] We found that at a larger interro- 
gation time (5ms) the fringe repetition rate becomes very 
high making accurate measurements of the fringe pat- 
tern difficult. Indeed, we found no clear pattern for the 
T = 5ms measurements, and the exact pattern for a given 
RF frequency was not perfectly repeatable. However, we 
note that qualitatively the distributions output by the 
interferometer at longer interrogation times are very dif- 
ferent from the results for T of order 100/is. Figs. [4jc) 
and (d) show column densities for a 5ms interrogation 
time in the case of a single tt/2 pulse and with the in- 
clusion of the second tt/2 pulse respectively. As Fig. [4^c) 
shows, any higher order processes do not make them- 
selves known in the distribution shape after 5ms, and 
indeed we have found that it takes almost an order of 
magnitude longer before mp changing collisions produce 
significant population in the mp = ±2 states. 

The result of applying a second tt/2 pulse after 5ms 
is shown in Fig. |4jd). The column density distribution 
in this typical case shows an essentially uniformly dis- 
tributed probability over the mp states. Such a distri- 
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bution is impossible to produce with coherent rotations 
alone. Dephasing due to the non-zero magnetic field gra- 
dient in the experiment is the likely cause of this behav- 
ior. We note in passing that the uniform distribution 
in Fig. [4^d) is reminiscent of spin transport results seen 
in [12] , where equipartition over spin sub-levels was found 
after a sufficiently long time. 

We now comment on the possibility of using different 
initial states as inputs to the interferometer. As seen 
in Figs. |4je-g), simulations show that merely changing 
the initially populated input port of the interferometer 
should produce non-trivial behavior in the fringes. In 
particular, taking |\E^) = |rap = +1) produces fringes 
in the mp = sublevel with four times the principle 
repetition rate of the interferometer. We have already 
measured the doubling of the repetition rate when the 
initial state is |rap = +2), but the usefulness of this ef- 
fect is debatable since the amplitude of the fringes is also 
halved in the mp = case. However, a quadrupling of 
the fringe repetition rate with no further loss in ampli- 
tude, as predicted here, should be useful since it provides 
four times the number of fringes to measure as compared 
to the usual case. 

In order to overcome the quantum limit of sensitivity, 
the current interferometer would in principle need to use 
squeezed input states. The relative ease of using tech- 
niques such as those of [5 to produce number-squeezed 
input states in our setup is promising, although calcula- 
tions are needed to confirm the precise benefits of such 
"twin atom" states for our interferometer. We would like 
to note, however, that even without considering exotic 
input states, improvements to sensitivity can be made, 
as seen in the predicted increase in the fringe repetition 
rate when the input state is in the = +1 magnetic 
sublevel (see the previous paragraph). Such increased ra- 
pidity of fringe oscillation compared with the fundamen- 
tal rotation frequency is often associated with entangled 
state interferometry [29] , but here arises naturally due to 
the multi-port nature of our interferometer. One motiva- 
tion for increasing the sensitivity is magnetometry: Be- 



cause /o depends on the local magnetic field, changes in 
B could be detected by looking for interferometric fringe 
shifts. 

We emphasize that the present results serve mainly 
to establish the possibility of high visibility Ramsey in- 
terferometry in a spin-2 Bose gas. Although the ob- 
served sensitivity is in line with results from other in- 
terferometry experiments, technical sources of noise need 
to be eliminated before the device can be used for pre- 
cision measurements. For interrogation times of order 
lms, dephasing effects come into play. In particular, 
the current experiment is limited by the non-zero mag- 
netic field gradient in the vacuum cell (~ 30mG/cm), 
and time-dependent fluctuations in the ambient magnetic 
field (amplitude ~ 100/iG). Such effects can be partially 
corrected for by using a spin echo 7r-pulse between the 
two Ramsey pulses used in the current experiment. In- 
deed, the overall goal for the present experiment is to sup- 
press sources of technical noise, and use such spin-echo 
pulse sequences for detection of AC magnetic fields [27] . 
As we will report elsewhere, such noise suppression and 
the application of more sophisticated pulse sequences can 
lengthen the interrogation time to several ms [28] . 

In conclusion, the atom interferometer realized here by 
applying coherent rotations to a spin-2 Bose gas, essen- 
tially functions as a multi-port interferometer with non- 
trivial fringes seen at each port. The measured fringe 
patterns demonstrate properties in good agreement with 
a phenomenological theory based on spinor rotation. The 
measured sensitivity is reasonable and may be improved 
firstly by eliminating technical noise in the imaging of 
the BEC, and ultimately by engineering the input state 
of the interferometer. Overall, our results provide a use- 
ful new tool for probing the spin coherent properties of 
Bose-Einstein condensates. 
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